Extended dynamical density functional theory for colloidal mixtures 

with temperature gradients 
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In the past decade, classical dynamical density functional theory (DDFT) has been developed 
and widely applied to the Brownian dynamics of interacting colloidal particles. One of the possible 
derivation routes of DDFT from the microscopic dynamics is via the Mori-Zwanzig-Forster projec- 
tion operator technique with slowly varying variables such as the one-particle density. Here, we use 
the projection operator approach to extend DDFT into various directions: first, we generalize DDFT 
toward mixtures of n different species of spherical colloidal particles. We show that there are in 
general nontrivial cross-coupling terms between the concentration fields and specify them explicitly 
for colloidal mixtures with pairwise hydrodynamic interactions. Secondly, we treat the energy den- 
sity as an additional slow variable and derive formal expressions for an extended DDFT containing 
also the energy density. The latter approach can in principle be applied to colloidal dynamics in a 
nonzero temperature gradient. For the case without hydrodynamic interactions the diffusion tensor 
is diagonal, while thermodiffusion - the dissipative cross-coupling term between energy density and 
concentration - is nonzero in this limit. With finite hydrodynamic interactions also cross-diffusion 
coefficients assume a finite value. We demonstrate that our results for the extended DDFT contain 
the transport coefficients in the hydrodynamic limit (long wavelengths, low frequencies) as a special 
case. 

PACS numbers: 82.70.Dd, 05.40.Jc, 05.45.-a, 47.57.E- 



I. INTRODUCTION 

While classical density functional theory has become 
a quite popular tool to calculate static properties of in- 
homogeneous fluids dHU, its generalization to dynam- 
ical, i.e., time-dependent, properties is much less ad- 
vanced!]] Most progress has been achieved for completely 
overdamped Brownian particles, which are realized as 
colloidal suspensions 0,11]. In a seminal paper of the 
year 1999, Marconi and Tarazona [H, [l(| have derived a 
dynamical density functional theory (DDFT) from the 
Langevin equations describing the motion of the individ- 
ual particles. The resulting DDFT equation corresponds 
to the field-theoretical model B for a single scalar order 
parameter, where the current is proportional to the func- 
tional density derivative of the equilibrium free-energy 
functional (generalized Fick's law). In 2004, Archer and 
Evans have used the stochastically equivalent Smolu- 
chowski picture to rederive this DDFT equation. In 
2009, Espahol and Lowen [l2| have employed the Mori- 
Zwanzig-Forster projection operator technique (MZFT) 
[l3l - [T9l | as a third derivation route by using the one- 
particle density as the only slow variable of the system. 

Subsequently, DDFT has been generalized toward bi- 
nary mixtures |20l - [23| and anisotropic particles (24l - [2r| 
as well as to the dynamics of freezing (2?i HI] and wet- 
ting [29i | . Moreover, solvent-mediated hydrodynamic in- 
teractions between colloids, which are typically neglected 



1 For a recent review, see Ref. [f|. 



in the modeling although they are important for actual 
colloidal samples, have been included into DDFT for the 
one-component case [30M32| . More recent generalizations 
concern particle self-diffusion in complex environments 
[33L HH , extern ally imposed flow fields [H, HI] , colloidal 
sedimentation [33], and "active" self-propelled particles 
[Mlll- 
In this paper, we follow the route via the MZFT in or- 
der to derive an extended DDFT (EDDFT), which goes 
beyond former DDFT in two respects. At first, we con- 
sider a multicomponent mixture of n different species of 
spherical, i.e., isotropic, colloidal particles. Hitherto, cal- 
culations for binary mixtures [2CH23| assumed a diagonal 
mobility matrix. Here, we show that there are in general 
nontrivial cross-coupling terms between the concentra- 
tion fields. We specify these non-diagonal terms explic- 
itly for colloidal mixtures with pairwise hydrodynamic 
interactions. Therefore, we establish the basic dynamical 
equations to apply DDFT to the dynamics of multicom- 
ponent colloidal systems including their hydrodynamic 
interactions. This constitutes a classic colloid problem, 
which has been explored intensely over severa l decades by 
using mode- coupling- like techniques 
simulations [H, Q , and experiments 



computer 
Secondly, 

we treat the energy density as an additional slow vari- 
able and derive formal expressions for an EDDFT con- 
taining also the energy density. This applies to situa- 
tions, where a nonzero temperature gradient is imposed 
leading to thermodiffusion, which is also known as the 
Ludwig-Soret effect (47l - l52j . The derived equations also 
incorporate the reciprocal effect, which is the so-called 
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Dufour effect |53| , where a concentration gradient causes 
energy transport. 

The MZFT is also the standard derivation route for 
mode-coupling theory (MCT) describing the dynamics 
of liquids. MCT has been applied both to molecular [54[ 
and colloidal liquids [55M58| and is by now an pivotal 
theory for the glass transition. Here, we are using the 
same theoretical framework of the MZFT to derive an 
EDDFT. Furthermore, we are using this common basis 
to compare EDDFT and MCT. A possible connection 
between DDFT and MCT has already been pointed out 
by Archer [HIIJ. 

There is also a close connection of the MZFT with 
several versions of the classical MCT close to phase tran- 
sitions [6ll - [64j |. which have been pioneered by Kawasaki. 
For example, for the isotropic-nematic phase transition 
in liquid crystal side-chain polymers a projection opera- 
tor technique of MZFT-type [65l [|36| has been used to 
construct a MCT [U HI|, which was able to explain 
the experimental results obtained by two groups [69147 1| 
as being due to a dynamic nonlinear coupling between 
order-parameter variations and the strain tensor. 

The use of MZFT to hydrodynamic condensed systems 
started with the work by Forster for nematic liquid crys- 
tals [HI, and has been applied in the following to a 
number of complex fluids with spontaneously broken con- 
tinuous symmetries 1 171 . in particular to the superfluid 
phases of 3 He I73l . l74j and, more recently, to uniaxial 
magnetic gels [75| . These applications were based on a 
generalization of the use of correlation functions in the 
hydrodynamic regime with applications to simple fluids 
[7q | and superfluid 4 He [72}. Thus, for all extensions 
of DDFT to more variables the hydrodynamic regime of 
long wavelengths and low frequencies emerges for all con- 
densed systems as a natural limit to check the results 
obtained. Conversely, the EDDFT can be used to inves- 
tigate how far the range of hydrodynamic considerations 
can be extended to larger values of frequencies and wave 
vectors. 

The paper is organized as follows: in Sec. [TTl we sum- 
marize in detail the technical aspects of the projection 
operator technique we use in a coherent fashion. In Sec. 
IIIIl we present the results of the application of the MZFT 
to colloidal mixtures in detail including the energy den- 
sity as a variable. Finally, we summarize our results 
and present a perspective for future generalizations of 
the present work in Sec. IIVI 



II. MORI-ZWANZIG-FORSTER TECHNIQUE 

The MZFT [U-[l5[ is described in detail in several text- 
books [l6l-[l~9l]. Further below, we comprehensively sum- 
marize the essential ideas that are relevant for this paper 
and adjust the notation to the problem at hand. 



A. General formalism 

For the purpose of this paper, it is most appropri- 
ate, but in general not necessary, to consider a grand- 
canonical ensemble of systems of iVparticles. The total 
ensemble T t with Hamiltonian i?(r t jj involves as canoni- 
cal variables the 6-/V coordinates qi(t) and momenta Pi(t) 
of the N particles. It can be described by the total proba- 
bility density p{t) = /5(r t ), which is given by the solution 
of the Liouville-von Neumann equation^ 

p = -tp = ~[H,p\, p(t) = e- £t m (1) 

with the Liouvillian £(T t ), the imaginary unit i, the re- 
duced Planck constant h = h/(2ir), and the commutator 
[A, Y] = XY — YX of A and Y. Alternatively, it is also 
possible to describe the same system in terms of only 
a few relevant variables a,; (t) = a,; (r, t) = a,; (f t ; fj3 with 
i = 1, . . . , n, which we assume to be real- valued in the fol- 
lowing. The corresponding relevant ensemble T t is asso- 
ciated with the relevant probability density p(t) = p(It). 
Using the relevant probability density p{t) and the grand- 
canonical trace Tr, which is given for classical systems by 

with the inverse thermal energy f3 = 1/ {k^T) 1 Boltzmann 
constant fee, absolute temperature T, chemical potential 
ji, and ensemble differential dT t = dqidpi • • ■ dq^wdpsN, 
the time-dependent ensemble average 

(A(0)) t = Tr(p(t)A(0)) = Tr(p(0)A(t)) (3) 

of an arbitrary variable X(t) can be defined. The aver- 
aged relevant variables a,i(t) = a,i(r, i jl are given by 

Oi(t) = (oj(0))t - Tr(ft(t)ai(0)) = Tr(ft(t)ai(p)) . (4) 



2 In order to keep the following expressions simple, the Hamil- 
tonian H(Tt) is assumed to be not explicitly time-dependent. 
Nevertheless, it is possible to consider also systems with a time- 
dependent external potential Ui(r,t) in the framework of the 
MZFT as it is described here, if the external potential varies suf- 
ficiently slowly with time so that it is approximately constant on 
microscopic time scales. 

3 In case of a classical system, the commutator [X, Y]/(ih) should 
be replaced by the Poisson brackets {X, Y} for any variables X 
and Y. 

4 The abbreviating notation Xi (t) = X; (r , t) = X; (Tt ; r) is 
used for several symbols X in this paper. Summation about 
an index i implies also integration over r for such symbols: 
XiYi = Jd 3 r Xi(r,t)Yi(r,t). (Einstein's sum convention is 
used throughout this paper.) 

5 Notice that a; (t) = ai (f, t) does not depend on the total ensemble 
ft, but the generalized summation rule defined in footnote [4] 
applies also here. 
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For a given thermodynamic functional like the Hclmholtz 
free-energy functional J 7 , their thermodynamic conju- 
gates a\(t) = a\{r,t) can be obtained by functional dif- 
ferentiation: 



5 2L 

Sat 



(5) 



A possible representation for the Helmholtz free-energy 
functional is 



1 



F = Tr(pH) + -Tr(pln(p)) 



(6) 



The Helmholtz free-energy functional .Ffajjl depends 
functionally on the averaged relevant variables cij(£) and 
is related to the grand-canonical functional ^[a* 1 ], that 
depends functionally on the thermodynamic conjugates 
a\ (t) , by the Legendre transformation 



b 

a] a; 



(7) 



A map from the total ensemble f t onto the relevant en- 
semble r t is constituted by a suitable projection operator 
Vt = 1 — Qt- This projection operator can be written as 

M 



P t X = Tr(p(t)X) + (a, 



(8) 



It projects onto a space that is spanned by the linearly 
independent variables di (t) and the unity id. 

The MZFT consists in the application of this operator 
in order to obtain transport equations for the relevant 
variables a,i(t), that are equivalent to the Liouville-von 
Neumann equations 



hi = Cat = - [H,ai] , 
n 



Oi(t) 



ai(0) , (9) 



by projecting out all irrelevant variables. These transport 
equations are given, for example, in Ref. [l6| in its general 
form. 

The dynamics of the reduced relevant variables 
Adi(t) = di(t) — a,i(t) is given by the (exact) general- 
ized Langevin equations 0, [I?} 



Aoj(f) = %(£)A%(i) + dt' Kijfat^Aajit') 

Jo (10) 



6 This functional JF[a\ = T[a\, . . . , a n ] is also called density func- 
tional in the context of DDFT. 

7 Equivalently, one could also choose a linearly dependent set of 
relevant variables a;(t) and construct a projector that maps onto 
the space spanned by these variables. In this case, the unity id 
is only indirectly taken into account as a basis element and the 
corresponding contribution - the first term on the right-hand- 
side of Eq. (|Sj - has to be omitted (see, for example, Ref. I Till ) . 



with the frequency matrix 



the memory matrix 







da 3 {t) 



Tr(p(t)a,i) 
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(12) 



and the noise 

Fi{t) = Q g(o,t)hi . (13) 

Here, Q(t',t) is the time-ordered exponential operator 



(14) 



where the time-ordering operator T- orders operators 
from left to right as time increases. 

Furthermore, the dynamics of the averaged relevant 
variables adt) is described by the averaged Langevin 
equations [16(] 



= Tr(p(t)hi) + fdt' Tv(p{t')C Q v g{t',t)hi) 
Jo 

+ F i (t) 



(15) 



with the averaged noise Fi(t) = Tr(p(0)i<i(i)). 

The frequency matrix fiy (£) takes the instantaneous 
reversible contributions to the dynamics of the relevant 
variables into account. In linearized form, it is an equal- 
time commutator of field operators. Whenever the cho- 
sen relevant variables di(t) have a definite time-reversal 
behavior, the frequency matrix vanishes. The memory 
matrix Kij(t,t'), on the other hand, comprises the non- 
instantaneous reversible contributions and all dissipative 
contributions to the dynamics of the relevant variables. 
The important finding that the memory matrix can also 
include (non-instantaneous) reversible contributions was 
first shown by Forster [l5|, [TtI [zH . 



B. Special generalized probability density 

The transport equations for the relevant variables and 
their correlation functions arc given in this section for 
the specific case of the generalized grand canonical prob- 
ability density 



Pit) = = 



1 



S(t) 



-l3H ott (t) 



(16) 



with the grand-canonical partition sum S(£) and the ef- 
fective Hamiltonian 



H eS (t) = H - 4(t)ai 



(17) 
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For this particular choice, the projection operator ((8j) is 
specified as the Robertson projectoid [HI, [78| 

P t X = Tr(p(t)X) + (a t -a J :(t))x^ 1 WTr(|^-x)(18) 
with the symmetric non-equilibrium susceptibility matrix 

Xij(t) = T~T7~: =0Tr(p(t)(ai - o< - a 3 -(f))) , 
Sa)(t) (19) 

the derivativ^l 

da\{t) 
and the operator 



= ^ t (a i -a j (*))p(t) , 



(20) 



(21) 



This operator can be omitted for a classical system: 
£ t X = X. 



1. N on- equilibrium dynamics 

The particular choice (|T5|) of p(t) leads to the exact 
transport equations [l6j 

di(t) = -By(t)o5(t) - Ut' Ri^t'^At') (22) 



with the antisymmetric drift matrix 

fl«(t) = l - H(p(t)[a j ,o J -]) = -B*(t) (23) 

and the retardation matrix 

%(i,t') =^Tr(p(t / )(4^(i , ,t)a i )(^a j )) . (24) 

Notice that these transport equations are applicable also 
far from thermodynamic equilibrium. 

2. Equilibrium correlations 

In thermodynamic equilibrium, the transport equa- 
tions for equilibrium time correlation functions (so-called 
Kubo functions [I?} ) 



C i3 {t) = (Aaf(*)|Aaf (0)> 



eq 



(25) 



with the equilibrium fluctuations Aa° q (t) = di(t) — a° q 
have no noise term. Here, the letters "eq" denote equi- 
librium quantities and Mori's scalar product is given by 



(X\Y) cq = Tr(p cq X<£ eq Y) 
with the equilibrium probability density 

1 



P 



eq 



and the operator 

£ eq Jf 



d\e-^ H Xe^ H 



(26) 



(27) 



(28) 



On the basis of Mori's scalar product (|26|) . the equilib- 
rium average (X) eq = {X\id) eq = Tr(p eq X) is defined. 
In the following, we present transport equations for the 
time correlation functions CV,(i). These transport equa- 
tions are at first given in position-time space (r , t) and 
later analyzed in the context of linear response theory in 
Fourier-Laplace space (fc, z). 

a. Position-time space: In the linear regime near 
equilibrium, the dynamics of the time correlation func- 
tions Cij (t) can be derived from Eq. (|22[) by linearization. 
The resulting transport equations are given by [HI, [l7| 



C ij {t) = Q^C kj {t) + dt'K%(t-t')C kj (t') (29) 
Jo 

with the equilibrium frequency matrix 



oeq eq-1 
n ik Afe 7 



and the equilibrium memory matrix 



(30) 



(31) 



These equilibrium matrices depend on the equilibrium 
drift matrisP] 



B if = ^([o» 5 Oj])eq 



-/3(Aar|A^ q ) 



on the equilibrium retardation matrbPl 

i%{t) = p{fr^(t)h i \h i ) eti , 

and on the static equilibrium susceptibility matrix 



(32) 



(33) 



(34) 



* For abbreviation, elements of inverse matrices are denoted as 

M" 1 = (M _1 )ij with an arbitrary matrix M in this paper. 
' This expression follows directly from the operator identity llfj . 

m 



da: 



Jo 



dAe AA W 



&A(x) 
Ax 



with a linear operator A(x). 



10 While the first equality in Eq. ((32} follows directly from Eq. 

the second equality holds only in the linear regime near equilib- 
rium and can be derived by linearization of the first term on the 
right-hand-side of Eq. H15H and comparison with the equivalent 
Eq. J22j. 

11 In Eq. I I 33 II . a redundant Q eq is often inserted directly in front 
of a; in order to obtain a more symmetric expression. 
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Here, the equilibrium projector V cq 
by 



1 — Q cq is given Here, the dynamic susceptibility matrix 



and the equilibrium exponential operator is 



g cq (t) 



(35) 



(36) 



Notice that the linearized transport equations (I29| can 
be used to determine the equilibrium frequency matrix 
fl° q and the equilibrium memory matrix K° q (t) exactly. 
It is thus possible to calculate these matrices by the eval- 
uation of equilibrium correlation functions obtained from 
experiments or microscopic simulations. 

b. Fourier- Laplace space: In Fourier-Laplace space 
(see appendix [XJ the dynamical equations (|2"9"|) obtain 
the simpler form [l7j 



\zd lk 



fl cq 



K 



cq 



A', 



k (z))C kj (z) = C ij (0) = J- 



with the obvious solution 



(37) 



(38) 



Notice that CV,-(z), C tj (0), fl® q , K° q (z) 7 and Xy are given 
in Fourier space, although their wave- vector dependence 
is not denoted explicitly here. Furthermore, X(t) denotes 
a time-dependent quantity, X(u>) its Fourier transform, 
and X{z) its Laplace transform in this paragraph. 

The Fourier transformed equilibrium frequency ma- 
trix fly and the Fourier-Laplace transformed equilibrium 

memory matrix K° q (z) in Eqs. (|37| and (|38| are given 
by 



fl oq 



with the Fourier transformed equilibrium drift matrix 



(40) 



the Fourier-Laplace transformed equilibrium retardation 



matrix 



lQ{z) = /3(A^ q |Q oc i(z - £^ q )- 1 Q cq |A^ q ) 



(41) 



with the equilibrium self-adjoined reduced Liouvillian 
= Q cq CQ cq [l^], and the Fourier transformed static 
equilibrium susceptibility matrix 



lim Yii(z)\ 



(42) 



~ r \ 1 (a *W 

71" JR UJ + 1Z 



and the absorptive response function 



X" j (t-t') = —(la t (t),a j (t')]) cq 



(43) 



(44) 



have been introduced. As usual in the context of lin- 
ear response theory, the absorptive response function ap- 
pears as a contribution in the complex response function 



(45) 



with the reactive part and the absorptive part 

X«(w), whose non-diagonal elements are not necessarily 
real- valued. The reactive response function xL(w) and 
the absorptive response function x'ij(u) are dependent 
and related to each other by the Kramers-Kronig (dis- 
persion) relations 



x:»^p/do/| 

K JR w 

x» : 



(46) 



Also the time correlation functions Cy(f) can be ex- 
pressed in terms of the absorptive response function: 



3y(*) = ^x£(*) • 



(47) 



Hence, their Laplace transforms Cij(z) are given by 



Cy (z) = 7^ / do; 



i-7r/3 J K w(w + iz) 



(48) 



C. Slow variables 

If the relevant variables vary sufficiently slowly with 
time so that there is a clear separation of time scales be- 
tween the slowly relaxing relevant variables and the fast 
relaxing irrelevant variables, the transport equations (|22[) 
and (|29l) can be simplified by neglecting contributions of 
order 0(af). Using the expansion pj| 



D Q t Ct = e Ct + ^ ^ 



(49) 



the retardation matrix (|24[) can be approximated by [T| 



12 The expression (14 1 I t has been symmetrized by insertion of redun- 
dant operators Q eq . 



, (t, t') = 13 Tr(p(t) (e^*-*') Q^i)^)) + O(Sjft50) 
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1. N on- equilibrium dynamics 



The approximation (|50[) results in the simplified trans- 
port equations for slow variables fl(| 

di(t) = -fl tf (t)o}(t) - /8i3y(t)o5(t) . (51) 

This Markovian approximation is also applicable far from 
thermodynamic equilibrium, but it is not appropriate, if 
effects related to "long time tails" (like the glass transi- 
tion [79|]) are investigated [ll|. The transport coefficients 
are given by the drift matrix Bij (t) and the mobility ma- 
trix Dij(t). The mobility matrix is given by the Green- 
Kubo-type expression 



D l3 (t) = Id? Ti(p(t)(e dt QA)(e 



(52) 



In the classical limit, this expression simplifies to fl2l.[l9| 

Dij(t) = Jd? Tr(p(t)(Q t a,)e £t '(4^)) , (53) 

where a redundant Qt has been inserted in front of dj in 
order to symmetrize the expression. Further redundant 
operators Qt could be inserted in the exponential func- 
tion in Eq. ([53)) by replacing the Liouvillian C by the 

self-adjoint reduced Liouvillian tf = QtCQt- Notice 
that the transport equations (|51[) in combination with 
the approximation (|53"|) are exact up to the third order 
in a,i. 



with local currents J^> (f, t) corresponding to at (r , t) . 
Analogous conservation laws hold for the averaged vari- 
ables Oj (r , t) with the averaged local currents (f, t) = 

Ti(p(t)J^(f,0)): a, + V ? - jW = 0. 

1. N on- equilibrium dynamics 

Since only classical systems with slow variables are 
considered in the following, dynamical equations for the 
time-evolution of the averaged relevant variables (r , t) 
can be derived from Eqs. ([IS]). (|5"T]) . ([53]). and ([ST]). These 
are the general classical extended DDFT equations 



cn(r,t) = - V r -Tr(p(t)J«(r,0)) 



n „ 

V V f • / dV PD<M (r ; r ', t)V P , a) (r \ t) 
, =1 Jm* 



(58) 



with the diffusion tensor 



di' Tt(p(t)(& J?\?, Q))e a \Q t J^{r, 0)) 



£t7 



f(0/ 



(59) 



If the variables dj(t) are real and have definite time- 
reversal signatures, one can show that D^p (r, r ', t) is 

symmetric 0: D { *p (r,r',t) = (rj^t). This state- 
ment is known as Onsager's principle [80|. 



,2. Equilibrium correlations 



2. Equilibrium correlations 

With the same approximation, the transport equations 
for the equilibrium time correlation functions Cij (t) 
become 

C ij (t) = 0£C kj (t) + r%C kj (fi). (54) 
Here, we introduced the transport matrix 

^=-mxir 1 (55) 

with the mobility matrix 

poo 

D% = d? (Aar|Q eq e £i 'Q eq |Aaf ) eq . (56) 



D. Conserved quantities 

An important example for slowly relaxing variables 
are local densities of conserved quantities. The trans- 
port equations of such conserved quantities hi (r , i) can 
be written as conservation laws 



The assumption of conserved quantities can also be 
used to rearrange the transport equations (|54[) for the 
equilibrium time correlation functions Cij (r , r ', t) into 

Cij (r, 7, t)=J2 Vr fd 3 r" L™ (r , f")C kj (r ", r\ t) (60) 
fc=i 3 
with the total transport matrix 

LM (r, r') = figf) (r, r') + (r, r') . (61) 

This matrix includes the contributions 



fci^ 3 (f 



(62) 



and 



a. + V^.jW =0 



(57) 



r^(r,r") = 2 d 3 r"PDW(r,r") V^f 1 ^') 
fe=i- yR3 (63) 

with the equilibrium drift tensor 

Bif (r, f') = P(J® (r, 0)| Afi^f, 0)) oq (64) 
and the equilibrium diffusion tensor 

poo 

D^{r,r') = /dt'(J^(r,0)|Q cq e £t Q ccl |J^(r",0)) cq . 

J ° (65) 
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III. 



COLLOIDAL MIXTURES 



We now consider a mixture of N c = Y17=i isotropic 
colloidal particles of n different species, where iV, is the 

total number of particles of species i G {1 rtjP^l 

These N c colloidal particles are suspended in a molec- 
ular solvent consisting of N small isotropic particles 
(molecules) of the same type. 

The MZFT is now used to derive an EDDFT equation 

for mixtures of colloidal particles. When r£' (t) denotes 

the position, p^\t) the momentum, and rrii the mass of 
the fcth particle of species i, where i = corresponds to 
the molecules of the molecular solvent and i > corre- 
sponds to the colloidal particles, the Hamiltonian of the 
system is given by 



N, 



#(r t ,i) = EE#n f *>*) 



(66) 



i=0 k=l 



with 



(67) 



3=0 1=1 



(r£\ t) is the external potential acting on the parti- 
cles of species i, U^if^ — fjf^) is the pair-interaction 
potential for two particles of species i and j, respectively, 
and r t is the total ensemble introduced in the beginning 
of Sec. [TTJ To assure that the MZFT as described in Sec. 
UT1 is applicable, the external potential is assumed to vary 
sufficiently slowly with time (see footnote [2J . The Liou- 
villian C(T t ,t) corresponding to the Hamiltonian (|66p of 
the considered system is 



£ = EE 



i=0 k=l 



(68) 



A. Relevant variables 

As relevant variables a,i(r,t) of the colloidal mixture, 
we choose the n concentrations 



,(f,i) = E^-4 i) (*)) 



(69) 



k=l 



with i G {1, . . . , n} and the energy density 

n Ni 



i=0 k=l 



13 Notice that the meaning of n in this section is different from its 
meaning in the previous section. 



Their averages are denoted as Ci(r,t) = Tr(p(0)ci(r, £)) 
and e(f,t) = Tr(p(0)e(r, t)) in the following. 

By considering only the n concentrations di (r , t) = 
Ci(r, t), . . . , a n (r, i) = c n {r,i) and the energy density 
a n +i(r, t) = i(r,t) as relevant variables, we assume that 

the momentum variables pi (t) relax much faster to local 
thermodynamic equilibrium than the position variables 
r^\t) so that the momentum density can be neglected 
as a further dynamic variable on the characteristic time 
scale of the concentrations and of the energy density. By 
this choice of relevant variables, we further assume that 
the concentration co(r, t) of the molecular solvent relaxes 
much faster than the concentrations Cj(f, t), i > 0, of the 
colloidal particles. 

The concentrations c, (f , t) and the energy density 
e(r, t) are even under parity and time reversal. Further- 
more, they are locally conserved, if there are no sources 
and sinks of particles and energy in the system. 

The corresponding currents follow from the Liouville 
equations a, + Iff, a,;} = [see Eq. by comparison 
with Eq. (|57pl 14 l They are the particle number current 



Ni Jj,) 



j c «(F,t) = 



(71) 



k=l 



and the energy current 



i=0 k=l % 

1 n Ni Nj Ji) Jj) 



i,j=0 k=l 2=1 



(72) 



xr^ d\5(r-rf + Arf>) 
Jo 



with the dyadic product <g> and the notation f 



kl 



ff , where all rjj. , r£j" , and H k in Eqs. (fTTj) and 
(1721) are to be taken at time t. 



Since J Ci (f 1 t) and J e (r, t) arc of odd order in the 
momentum p~ k l \t), the averages Tr(p(£)J Ci ) = and 
Tr(p(i)J e ) = vanish. This leads to the important in- 
variancc properties QtJ Ci = J Ci and QtJ £ — J e [see Eq. 



For the derivation of J E (r , t ) , the equation [ltj ] 
S(f~r')-S(f-r") 

= -V r ((?' - r") JdX <5(r - f" + A(r ' - f "))) 

is helpful. 
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B. Transport equations 



These dynamical equations are given by 



Since the concentration fields Ci (r , t) and the energy 
density e(r, t) are locally conserved, the EDDFT equa- 
tions (|58|) and the corresponding transport equations 
respectively, can be applied. Due to the invari- 
ance of c\(r,t) and i{r,t) under time-reversal, the fre- 
quency matrix and therefore also the first term on the 
right-hand-side of Eq. (|58|) vanish. 



1. N on- equilibrium dynamics 

Application of Eqs. ([55)1 results in the following ex- 
tended DDFT equations for colloidal mixtures: 

n „ 

ti (r ; t) = Vv f d V (r, r ', t) V f / c5 (r ', t) 

^ (73) 



+ V r ~ / d 3 / pD<& (r, f', t)V r -,£ h (r', t) , 

n „ 

U ^ (74 ) 
+ V r - / dV f3D^ (r, r", ijVjwe^f , *) . 



The diffusion tensors £)W) (r*, f', i), D^(r,f',t) = 
(D < - £ ^(r,r',t)) T , and (f,f',t) in the EDDFT equa- 
tions (|75|) and (f?4")) are given by 

/■CXD 

D™ = jdt> Tr(p(t)J^(r,t')J^(f',0)) , (75) 

/■CXD 

(r, r", t) = Jdt' Tr(p(t)4- (r, t')Jf{f', 0)) , (76) 

/■CXD 

Dj^ = /df TVUt)J|(r,tOJf(r',0)) . (77) 



They are associated with particle diffusion {D^>), (in- 
verse) thermodiffusion (D^: Ludwig-Soret effect, D^: 
Dufour effect), and heat conduction (Z)( e£ )), respectively. 
The EDDFT equations ([73]) and (J7J]) in combination with 
the diffusion tensors (j75|) - (f7Tj) constitute the main result 
of this paper. 



2. Equilibrium correlations 

If the transport equations (|60[) arc applied, one obtains 
dynamical equations for the time correlation functions 

GjftCt) = (Ac?(r,t)\Ac?(r\0)) cq , (78) 
C ie (r, f', t) = (AcJV, t)\Ai cq (r\ 0)) cq , (79) 
C eE (r, r\ t) = (Ae^(r, t)\Ae^(r', 0)} cq . (80) 



fc=i 



Ga (r- , r ', t) = ^V r - d V L W (r, r ") C kJ (r ", r ', t) 



+ V r /d 3 r" (f, r")C ej (r ", f ', t) 



(81) 



n „ 

C J£ (r , r ', i) = V V r / d 3 r" L &*> (r, r ") C fe£ (r ", r ', t) 
fe=i 7r3 



+ V ? - /dV' 4-)(f, r")C ££ (r", r', i) 



(82) 



n » 

C e£ (r , r ', i) = V V r / d 3 r" L&*> (r , r ") C fe£ (r ", r ', t) 
fc=i 7r3 



+ V ? - /dV' 4f(f, r")C ££ (r", f', i) 



(83) 



with the total transport matrices 

L«> (r, f ') = fig) (r, r") + rW) (r, r') , (84) 

^' ?) = *") + It £) (r, r") , (85) 

Lif (r, r') = nif (r, f) + lf q e ) (r, r') (86) 

consisting of the contributions 



fc=i 



d 3 r"B(-)(r,f")x:r 1 (f",f'), 



(87) 



n&'c.o =-E ^"BM^nxtrH^f') 

fe=l R3 



dV'M^^r'OxS- 1 ^''^) 



(88) 



n /. 
h=i • /Ra 



(89) 



R 3 



and 



n « 

-^ R3 (90) 

(9i 

+ /dV^^^OV^xS- 1 ^''^') , 

JR3 

n « 

k=i J ^ 3 



(91) 



(92) 



+ / d 3 r" f3Dif (f, r ") V^xS"^^' r') 



R3 
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Equations (|87 |) -(|92 p in turn depend on the equilibrium 
drift tensors 

B^\r,f') = /?(>^,0)|Acf (f',0)> eq , (93) 

(r, f) = p(h (r, 0)\A^{f', 0)} cq , (94) 

B ( c f (r, r') = /?{> (r, 0) | Ae^(r ', 0)) eq , (95) 
on the equilibrium diffusion tensors 

D^(f,n = dt' (J c '(r,0)|^(r',t'))e q , (96) 
Jo 

(r, 7) = / dt' (J c ' (r, 0) | J' (f ', t')) cq , (97) 
Jo 

/■OO 

(98) 



Dit e) (f,n = dt' (J%f,0)\J £ (r',t% 



and on the static equilibrium susceptibility matrices 

X-f(f, r') = (3(A^(r, 0)|Acf (f , 0)) eq , (99) 
X ^(r, r') = (3(A£°\r, 0)| A^(f ' 0)} cq , (100) 
X° q (r,r") = /3(A£ oq (r,0)|Ae cq (r',0)) oq (101) 

with the equilibrium fluctuations Ac° q (r, t) = Ci(r,t) — 
c| q (r) and Ae cc i(r, t) = e(r, t) - e e %r). 

C. Approximation of the diffusion tensors 



For an application of the EDDFT equations ([73]) and 
(|74[) to a particular system, suitable expressions for the 
diffusion tensors (fT5"j) - ([7"T|) are needed. A possibility to 
determine these diffusion tensors is the implementation 
of particle- resolved computer simulations |8l| . Alterna- 
tively, analytical approximations for the diffusion tensors 
D^(r,r',t), D^ ie \r,r',t), and D^ e \f,r',t) can be ap- 
plied. Such approximate expressions are given in the fol- 
lowing. 



1. No hydrodynamic interactions 

As first approximation, it is assumed that the consid- 
ered system is sufficiently close to local thermodynamic 
equilibrium so that the relevant probability density p(t) 
can be approximated by the equilibrium probability den- 
sity p cq [see Eq. ([27])] in Eqs. (|75 ]) -(|77 ) . Secondly, we as- 
sume that the position variables relax much more slowly 
to local thermodynamic equilibrium than the momen- 
tum variables and that the external potential is approx- 
imately constant on microscopic length scales. Thirdly, 
we suppose that the position and momentum variables 
are statistically independent. Furthermore, the consid- 
ered suspension shall be sufficiently dilute so that hydro- 
dynamic interactions between the colloidal particles can 
be neglected and the momenta of different particles are 
uncorrelatcd. Finally, we assume oricntational isotropy 



for the momentum variables, i.e., {p®p)eq = ^l(P'P)cq 

and neglect the pair-interaction potential \^ —fj? ) 
in Eqs. (|S7|) and 

With these assumptions, the diffusion tensors (|75p - ([77|) 
can be approximated by 

D™(?,r',t) = D$hd ij 6(?-f%(r,t) , (102) 
D&l (r, r', t) = Kf 16(r- r') Cl (r , t) , (103) 

n 

^nh (r, f', t)=J2 «H lS(f- f)a (r, t) (I 04) 

i=0 



with the transport coefficients 

D^ = -\dt' (v i (t')-^(0)) eq , 



.(') 



3 Jo 



dt' (v^f) -^(0)^(0)), 



1 c°° 

= -Jdf (VitfyffiWHiWHM)^ 



(105) 
(106) 
(107) 



where Vi (t) = p^^ (t) / rrii is the velocity of a colloidal par- 
ticle of species i and Hi(t) is its energy. These coeffi- 
cients are associated with particle diffusion, thermodiffu- 
sion, and heat conduction, respectively Notice that the 
diffusion tensor (|102l) is diagonal and that all diffusive 
cross-couplings in Eq. (|73|) vanish, if there arc no hydro- 
dynamic interactions between the colloidal particles. 



2. Hydrodynamic interactions 

A better approximation for the diffusion tensor (|75|) . 
that takes also diffusive cross-couplings into account, can 
be derived, if hydrodynamic interactions between the 
colloidal particles are taken into account. In order to 
do so, the derivation of the DDFT equation for a one- 
component suspension of colloidal particles with hydro- 
dynamic interactions in Refs. [3(| El is generalized and 
compared with Eq. ([73)) . This derivation starts from the 
Smoluchowski equation [l8[ 

n Ni 

p(? N > f )+EE 4:l(^» *) = o (108) 

1=1 k=l k 



with the iV-particle probability density P(r N ,t), where 

^ = (rf ) ,...,r-g) 1 
probability currents 



f N = (fi , . . . , ) are the positions of all particles, the 



n N 3 

4%(r N ,t) = -EE D « , (109) 

j=l 1=1 

and the force densities 

(f N , t) = P(r N , t) + P(r N , i)V^ (f3U(r N , t)) . 



(110) 
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Here, Dj^r*) is a short-time diffusion tensor and 
U(r N 7 t) denotes the totai potentiai energy of the sys- 
tem. If the considered suspension is not too dense so 
that the particle distances are sufficiently large, the hy- 
drodynamic interactions can be approximated on the 
two-particle level and higher-order hydrodynamic inter- 
actions are negligible. 

In case of only hydrodynamic pair-interactions, the 
short-time diffusion tensors T)^p{f N ) can be written in 
the exact form [H H3| 

DjW(i*) = 0$V*l 

n N q 

+ d^Mh - S M hi s iq) (*t 9) ) ( m ) 

3=1 p=l 

+ (i-%fe)^' ) h^ ) (rif) 

with the self- and cross-interaction functions 

h ( * j) (r) = A^ j) (r)r®f + J B^' ) (r)(l - f®f) (112) 

with A = s for "self and A = c for "cross" , respectively, 
and the notation r = ||F|| and f = f/r for an arbitrary 
vector r. The self- and cross-interaction functions de- 
pend on the four mobility functions (r) and b[ 1 ^ (r) 
with A G {s,c}. With the method of reflections [18l l82j. 
these mobility functions can be determined up to arbi- 
trary order as an expansion in the inverse inter-particle 
distances. Up to fourth order, the mobility functions are 
given bJ^l 

A«)(r) = C9(r- 4 ), (113) 
B p( r ) = o(r- 4 ), (114) 

B^r) = \^ + \^^ + G(r-±), dl6) 

where R t denotes the radius of a colloidal particle of 
species i. Notice that Eqs. (|111[) and f)l 1 2[) together with 
the fourth-order approximations ()113[) - (|116[) of the mo- 
bility functions constitute ^generalized Rotne-Prager ap- 
proximation for mixtures flo . |83| . 

The generalization of the derivation in Refs. [3(J EJ 
leads to the following approximation of the diffusion ten- 
sor (|75[) for hydrodynamic pair-interactions: 

D^(f,r',t) = DjpSijSif- n^Ci(r,t)+c^(f,t)) 

+ D^\r-?')c ij {r,r',t) . ( 117 ) 



Here, we introduced the functions 

n « 

<4V>*) = V /dV h»)(f- r') Cij {r,f",t) (118) 

and the two-particle concentrations 

Cij (r, r', t) = Tr(p(0)c„- (f, r ', i)) (119) 
with the corresponding variables 

Ni Nj 

(r, r ', i) = EE^- *f (*)) ^ ' - ^ (*)) ■ ( 12 °) 
fc=l !=1 
l^k 

The two-particle variables (r , r ', t) are assumed to re- 
lax much faster to local thermodynamic equilibrium than 
Ci(r, t) and e(r, t). 



D. Approximation of the free-energy functional 

In order to determine the Hclmholtz free-energy func- 
tional J-[ci, . . . , c n , s], which is needed in the EDDFT 
equations ([73]) and (|74|) . it is always possible to expand 
this functional with respect to Ci(r, t), e(r,t), and their 
gradients taking general symmetry considerations into 
account [Hj|. 

If the energy density can be neglected so that only an 
approximation for the functional J-[ci, . . . , c n ] is needed, 
static density functional theory can be applied to derive 
such an approximation on a microscopic basis. The up to 
now most accurate approximation for J-[c\ , . . . , c n ] was 
derived in the framework of fundamental measure theory 
(see Ref. [H| for a review). 



E. Special cases of the EDDFT equations 

The EDDFT equations ([73]) and ([74]) contain several 
special cases that are known from the literature or that 
are relevant for particular applications. Two of these spe- 
cial cases are addressed in this section. The first one is 
an isothermal binary mixture, where only two concen- 
trations are present and the energy density can be ne- 
glected. As a second example, the hydrodynamic limit of 
the EDDFT equations is discussed. 



1. Isothermal binary mixture 



« The identity A^>(f£f ®r£f) = (21- 6 f <f ® r«>)/r«> 2 is 

' k 

useful for the derivation of Eqs. H113t -l lll6ll . 



If the considered mixture consists only of n = 2 differ- 
ent species of colloidal particles and the energy density 
can be assumed to be constant, the EDDFT equations 
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(121) 



(|73|) and (|74|) can be simplified to 

c 1 (f , = W V f • (1 Cl (r , t) + c« (r , t)) V f (r , t) 
+ PD$>Vr / d 3 r' h( n )(r - f') Cll (r, f, t)V f ,c\(r', t) 
+ pD$> Vr- 1 dV h( 12 > (r -f')c 12 (r, r', t)V F ,c\ (r', t) 

JR 3 

c 2 (r,t) = ^ 2) V r -(lc 2 (r,t) +c^(r,t))V?c\{r 1 t) 

+ pD^Vp- [ dV h^\f- r')c 21 (r, r', t)V F ,c\(r', t) 
Jr 3 

+ PD^V F - f dV h< 22 ' (r - f ')c 22 (r, f', t)V^c 2 (f', t) 



with the constant diffusion coefficients 



(122) 



Here, the short-time diffusion coefficients (|105[) can be 
expressed by 



D 



1 



(123) 



with the dynamic (shear) viscosity r\ of the molecular sol- 
vent. Hydrodynamic interactions between the colloidal 
particles are still taken into account by Eqs. (jT3")l and 
(|74|) . In the fourth-order approximation (|113p - (|116|) . the 



functions ci (r, t) = vanish (0 denotes the zero matrix) 



and the cross-interaction functions hi (r) are 



4||r|| V IMIV 



1 i? 2 i?j + R 3 



1-3 



(124) 



As closure relations for the two-particle concentrations 
Cij(r,r',t) in the dynamical equations (|121j) . the (ex- 
act) generalized Ornstcin-Zcrnikc equation for mixtures 
or simple analytical approximations that are known from 
the literature can be applied [3(1 HH . 



2. The hydrodynamic limit 



The derived EDDFT equations (J73J) and ([74]) with 
the space- and time-dependent diffusion tensors (|75l) - ([77|) 
constitute an extension of the corresponding hydrody- 
namic equations to larger wave vectors k and frequencies 
ui. In the hydrodynamic limit (k —> 0, ui — > 0), the 
EDDFT equations become 

n 

(f , t) =£ )A ? 4 ^ *) + P°0 ^ ft *) . 

3=1 (125) 

n 

e(f, i) = £ /3^' } A r -4^ *) + 

3=1 (126) 



W) 



d 3 r At Tr(p(0)J Cl (f,t)-J c '(0,0) , (127) 



D { ee) = i /d 3 r 



p(0)>(r,t)->(d,0)) (129) 



and Dq"'' = D^"\ The hydrodynamic limit of the trans- 
port equations (|8T |) - ((83|) for the time correlation func- 
tions ((75)) - ([50)1 can be obtained analogously. 



F. Relation of EDDFT and MCT 

The MCT of glass transitions [54| [79| is a classical the- 
ory for the dynamics of liquids near the glass transition. 
Originally, MCT was constructed for the underdamped 
dynamics of atomic and molecular systems [79j], but it 
can also be derived for the overdamped dynamics of col- 
loidal systems [H, H3, HI] . Like EDDFT, also MCT can 
be derived from the MZFT. This allows a comparison of 
these two theories on a common fundamental basis. In 
the following, we summarize the derivation of MCT and 
discuss its relation to DDFT and EDDFT. 



1. MCT for atomic and molecular systems 

The traditional form of MCT applies to a one- 
component system of equal spherical atoms or molecules 
of mass to. This system is characterized by a one- 
particle density field c(r , t) following the conservation law 

c(r, t) + V?- J c (r, t) = with the density current J c (r, t). 
In order to derive MCT, we switch to the Fourier-La- 
place space and utilize Eqs. (|37p . where we omit the let- 
ters "eq" denoting equilibrium quantities and the tilde ~ 
denoting quantities in the Fourier-Laplace space for rea- 
sons of clarity in this section. Near the glass transition, 
two variables are taken into account as relevant variables. 
These are the density field c(fc, z) and the longitudinal 

component j L (fc, z) = k/k ■ J c (k,z) of the density cur- 
rent J c (k,z). The transversal component of the density 
current, on the other hand, does not couple to density 
fluctuations and can therefore be neglected. While this 
was not the case in the context of EDDFT, here also 
the current associated with the density field has to be 
regarded as a relevant variable, since there is no sepa- 
ration of time scales between these variables near the 
glass transition [79[ . We further define the concentration 
time autocorrelation function (dynamic structure factor) 
C c (k, z) = (Ac(k, z)\Ac(k, 0)} and the current time auto- 
correlation function C 3 (k, z) = (Aj L (k, z)\Aj L (k, 0)) cor- 
responding to the chosen relevant variables. With these 
definitions, application of Eqs. (|3"T[) leads directly to a 
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dynamical equation for the normalized density time au- 
tocorrelation function 4> c (k,z) = C c (k, z)/C c (k, 0). This 
dynamical equation is the MCT equation [79| 

-l 



c (k,z) =[z + 



Ki{k,z) 



(130) 



with the frequency f2 m (/c) that must not be confused 
with the (vanishing) frequency matrix in Eqs. (|37j) . This 
frequency is given by fi^(fc) = C j (k, 0)fc 2 /C c (fc, 0) and 
C J (fc, 0) = 1/(/3to). Furthermore, the current memory 
function K^(k,z) in Eq. (|130[) is defined as K^(k,z) = 

-(A^(k,0)\Q(z - Ctf-iQl^fcoy/CPfco). 

2. MCT for colloidal systems 

In case of a system of spherical colloidal particles that 
are suspended in a molecular solvent, a simpler MCT 
equation but with the same long-time behavior as Eq. 
(|130p can be derived. This colloidal system is charac- 
terized by the short-time diffusion coefficient Do and 
concentration field c(k, z) of the colloidal particles. A 
similar derivation as before, but now with the appro- 
priate Smoluchowski operator £§ instead of the Liouvil- 
lian £, leads to the MCT equation for colloidal systems 



4> c (k,z) = [z + 



nl(k) 



1 - Ki(k,z) 



(131) 



with n*(fc) = D k 2 /C c {k,0). Notice that Eq. (fTBTjl is 
only of first order in z, while Eq. (|130p is of second order. 



3. Comparison of EDDFT and MCT 

Although the derivation of MCT was only presented 
for the simple special case of a one-component system 
here, more general formulations of MCT exist that are 
like EDDFT, for example, also applicable to (colloidal) 
mixtures [79j . Even the incorporation of the energy den- 
sity into MCT has already been discussed [86| in the lit- 
erature. EDDFT and MCT are therefore two different 
general theories with overlapping fields of application. A 
possible relation of DDFT and MCT has been mentioned 
by Archer [5!^, Ho[, but was not yet rigorously proven. 
Archer showed that under certain approximations the 
traditional DDFT equation f9j-fll| can be rearranged into 
a transport equation for the density time autocorrelation 
function, which matches the standard form (|131l) of MCT 
for colloidal systems. However, his derivation, which sug- 
gests that MCT can be derived from DDFT, is not rigor- 
ous, since it involves a reinterpretation of the one-particle 
density field as a temporally coarse-grained density field. 

In contrast, our derivation of the EDDFT presented in 
this paper allows to compare both theories from a fun- 
damental point of view. The derivation of EDDFT and 



MCT on the basis of the MZFT is illustrated in Fig. [Q 
A comparison of the derivations of these theories makes 
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FIG. 1. Illustration of the derivation of EDDFT and MCT 
using the MZFT. 

clear that there are actually strong differences between 
EDDFT and MCT and that a rigorous derivation of one 
theory from the other is not possible. An obvious and 
important difference between EDDFT and MCT results 
from the different approximations made in their deriva- 
tions. While the derivation of EDDFT involves a Marko- 
vian approximation when slow variables are assumed (see 
Sec. Ill C 1[) . the MCT equations are non-Markovian - a 
feature that becomes indeed relevant near the glass tran- 
sition. A more detailed comparison reveals that EDDFT 
and MCT are rather complementary but not replaceable 
theories. While EDDFT has proven that it is successfully 
applicable to describe weakly correlated systems with low 
concentrations, it cannot be applied to the glassy dynam- 
ics of systems at very high densities, since the Markovian 
approximation in the derivation of the EDDFT equation 
can only be justified, if effects associated with long time 
tails can be neglected [l6|. To the contrary, the deriva- 
tion of MCT does not involve a Markovian approxima- 
tion and has proven to be a useful analytical tool for the 
description of strongly correlated systems with high con- 
centrations, where long time tails have to be taken into 
account [79(- However, its derivation involves strong ap- 
proximations, too, so that MCT fails when it is applied 
to weakly correlated dilute suspensions. 



IV. CONCLUSIONS 

In this paper we have generalized classical dynami- 
cal density functional theory (DDFT) using the Mori- 
Zwanzig-Forster projection operator technique (MZFT) 
by adding concentration fields and the energy density 
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as variables. The resulting extended dynamical density 
functional theory (EDDFT) was compared to its hydro- 
dynamic limit and to mode-coupling theory (MCT) re- 
vealing that EDDFT and MCT are complementary the- 
ories with different fields of application. Our EDDFT 
framework shows that the MZFT is a flexible framework 
to incorporate thermal gradients (and other possible slow 
fields). 

We emphasize that, in principle, our EDDFT equa- 
tions ([75)1 and ([71]) treat concentration and temperature 
gradients on arbitrary length scales even down to micro- 
scopic length scales of the average distance between the 
colloidal particles. The essential input for our EDDFT 
equations are functional derivatives, which can be ob- 
tained from equilibrium correlations, and diffusion ten- 
sors, which can be obtained from dynamical correlations. 
An important challenge for the future is to apply this 
concept to actual temperature gradients in order to pre- 
dict the Soret coefficient. 

Guided by the application of the MZFT to various 
hydrodynamic systems including those with macroscopic 
degrees of freedom associated with spontaneously broken 
continuous symmetries [HI, [zl, [74| , it will also be inter- 
esting to see to what extent one can generalize hydrody- 
namic considerations [76l |77| using correlation functions 
to larger wave vectors and frequencies. 

A future generalization of the EDDFT equations 
should also take anisotropic colloidal particles with 
macroscopic degrees of freedom into account so that col- 
loidal liquid crystals can be addressed [2(|. It will be 
important to compare such an approach to the results 
obtained previously for colloidal liquid crystals using a 
paramctrization of the density with spherically symmet- 
ric, dipolar, and quadrupolar contributions [ 



Also the incorporation of the entropy density [9CJJ as a 
further variable would be an important task for the fu- 
ture. 

Another challenge for the future is the potential use of 
the MZFT for systems driven far from thermodynamic 
equilibrium for which a generalized thermodynamic po- 
tential }9ll - l94| is not known. To address this question 
appears to be particularly important for active systems, 
which have increasingly come into focus over the last few 
years ^MOk- 

Recently, a similar approach using the MZFT for the 
one-particle densit y an d the energy density was put for- 
ward by Espahol [10C| . This approach is based on an 
entropy functional formalism and provides explicit ex- 
pressions for hard spheres. However, mixtures and hy- 
drodynam ic in teractions are not treated explicitly in this 
approach [l00l |. 
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Appendix A: Integral transformations 

Since there are different definitions of the Fourier- and 
Laplace transformations in the literature, here we sum- 
marize the definitions that have been used in the con- 
text of the work presented. In addition, two useful re- 
lations between the Fourier- and Laplace transformation 
are given. 

1. Fourier transformation 

The Fourier transformation of a space- and time- 
dependent function X(r,t) is given by 



X(k,uj)= /d 3 r j&tXir^y^-^ 

1 



X(r,t) 



(2tt)< 



d 3 k dwll^le" 1 ^ 



(Al) 



with k € R 3 and 



2. Laplace transformation 

The Laplace transformation of a time-dependent func- 
tion X(t) is given by 



X(z) = dtX(t)e- zt , 
Jo 



1 



c+ioo 



(A2) 



X(t) = — ldzX(z)e zt 

27rl Jc_ioo 



with z£C and the real part > 0. The expression for 
the inverse Laplace transformation is known as Bromwich 
integral and contains a constant c > zq, where zq is the 
convergence abscissa of X(z). 

3. Useful relations 

The Fourier transformed function X(lo) and the 
Laplace transformed function X(z) can directly be trans- 
formed into each other. With the residue theorem, the 
following map from X(u) to X(z) can be proven: 



2nJ m to + iz 



(A3) 



A complementary map from X(z) to X(oS) is given by 

(A4) 



:(u>) = lim (x(z)\ -X(z) 



z—iu — e 



It follows di rectl y from the (special) Sokhotski-Plemelj 
theorem [l7l.[Toi| 



(A5) 



lim — = P— ± inS(x) , 

e->o+ x =F ie x 

where P denotes the Cauchy principal value. 
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